An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

The secreted triose phosphate isomerase of Brugia malayi is required to sustain microfilaria production in vivo

Human lymphatic filariasis is a major tropical disease transmitted through mosquito vectors which take up microfilarial larvae from the blood of infected subjects. Microfilariae are produced by long-lived adult parasites, which also release a suite of excretory-secretory products that have recently been subject to in-depth proteomic analysis. Surprisingly, the most abundant secreted protein of adult Brugia malayi is triose phosphate isomerase (TPI), a glycolytic enzyme usually associated with the cytosol. We now show that while TPI is a prominent target of the antibody response to infection, there is little antibody-mediated inhibition of catalytic activity by polyclonal sera. We generated a panel of twenty-three anti-TPI monoclonal antibodies and found only two were able to block TPI enzymatic activity. Immunisation of jirds with B. malayi TPI, or mice with the homologous protein from the rodent filaria Litomosoides sigmodontis, failed to induce neutralising antibodies or protective immunity. In contrast, passive transfer of neutralising monoclonal antibody to mice prior to implantation with adult B. malayi resulted in 60–70% reductions in microfilarial levels in vivo and both oocyte and microfilarial production by individual adult females. The loss of fecundity was accompanied by reduced IFNγ expression by CD4+ T cells and a higher proportion of macrophages at the site of infection. Thus, enzymatically active TPI plays an important role in the transmission cycle of B. malayi filarial parasites and is identified as a potential target for immunological and pharmacological intervention against filarial infections

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

Molecular dynamics simulations of the growth of poly(chloro-para-xylylene) films

Parylene C, poly(chloro-para-xylylene) is the most widely used member of the parylene family due to its excellent chemical and physical properties. In this work we analyzed the formation of the parylene C film using molecular mechanics and molecular dynamics methods. A five unit chain is necessary to create a stable hydrophobic cluster and to adhere to a covered surface. Two scenarios were deemed to take place. The obtained results are consistent with a polymer film scaling growth mechanism and contribute to the description of the dynamic growth of the parylene C polymer

Preference of Small Molecules for Local Minimum Conformations when Binding to Proteins

It is well known that small molecules (ligands) do not necessarily adopt their lowest potential energy conformations when binding to proteins. Analyses of protein-bound ligand crystal structures have reportedly shown that many of them do not even adopt the conformations at local minima of their potential energy surfaces (local minimum conformations). The results of these analyses raise a concern regarding the validity of virtual screening methods that use ligands in local minimum conformations. Here we report a normal-mode-analysis (NMA) study of 100 crystal structures of protein-bound ligands. Our data show that the energy minimization of a ligand alone does not automatically stop at a local minimum conformation if the minimum of the potential energy surface is shallow, thus leading to the folding of the ligand. Furthermore, our data show that all 100 ligand conformations in their protein-bound ligand crystal structures are nearly identical to their local minimum conformations obtained from NMA-monitored energy minimization, suggesting that ligands prefer to adopt local minimum conformations when binding to proteins. These results both support virtual screening methods that use ligands in local minimum conformations and caution about possible adverse effect of excessive energy minimization when generating a database of ligand conformations for virtual screening

Normal-Mode-Analysis–Monitored Energy Minimization Procedure for Generating Small–Molecule Bound Conformations

The energy minimization of a small molecule alone does not automatically stop at a local minimum of the potential energy surface of the molecule if the minimum is shallow, thus leading to folding of the molecule and consequently hampering the generation of the bound conformation of a guest in the absence of its host. This questions the practicality of virtual screening methods that use conformations at local minima of their potential energy surfaces (local minimum conformations) as potential bound conformations. Here we report a normal-mode-analysis–monitored energy minimization (NEM) procedure that generates local minimum conformations as potential bound conformations. Of 22 selected guest–host complex crystal structures with guest structures possessing up to four rotatable bonds, all complexes were reproduced, with guest mass–weighted root mean square deviations of <1.0 Å, through docking with the NEM–generated guest local minimum conformations. An analysis of the potential energies of these local minimum conformations showed that 22 (100%), 18 (82%), 16 (73%), and 12 (55%) of the 22 guest bound conformations in the crystal structures had conformational strain energies of less than or equal to 3.8, 2.0, 0.6, and 0.0 kcal/mol, respectively. These results suggest that (1) the NEM procedure can generate small–molecule bound conformations, and (2) guests adopt low-strain–energy conformations for complexation, thus supporting the virtual screening methods that use local minimum conformations

The impact of helminths on the response to immunization and on the incidence of infection and disease in childhood in Uganda: design of a randomized, double-blind, placebo-controlled, factorial trial of deworming interventions delivered in pregnancy and early childhood [ISRCTN32849447]

BACKGROUND: Helminths have profound effects on the immune response, allowing long-term survival of parasites with minimal damage to the host. Some of these effects "spill-over", altering responses to non-helminth antigens or allergens. It is suggested that this may lead to impaired responses to immunizations and infections, while conferring benefits against inflammatory responses in allergic and autoimmune disease. These effects might develop in utero, through exposure to maternal helminth infections, or through direct exposure in later life. PURPOSE: To determine the effects of helminths and their treatment in pregnancy and in young children on immunological and disease outcomes in childhood. METHODS: The trial has three randomized, double-blind, placebo-controlled interventions at two times, in two people: a pregnant woman and her child. Pregnant women are randomized to albendazole or placebo and praziquantel or placebo. At age 15 months their children are randomized to three-monthly albendazole or placebo, to continue to age five years. The proposed designation for this sequence of interventions is a 2 x 2(x2) factorial design. Children are immunized with BCG and against polio, Diphtheria, tetanus, Pertussis, Haemophilus, hepatitis B and measles. Primary immunological outcomes are responses to BCG antigens and tetanus toxoid in whole blood cytokine assays and antibody assays at one, three and five years of age. Primary disease outcomes are incidence of malaria, pneumonia, diarrhoea, tuberculosis, measles, vertical HIV transmission, and atopic disease episodes, measured at clinic visits and twice-monthly home visits. Effects on anaemia, growth and intellectual development are also assessed. CONCLUSION: This trial, with a novel design comprising related interventions in pregnant women and their offspring, is the first to examine effects of helminths and their treatment in pregnancy and early childhood on immunological, infectious disease and allergic disease outcomes. The results will enhance understanding of both detrimental and beneficial effects of helminth infection and inform policy

Conformational Control of the Binding of the Transactivation Domain of the MLL Protein and c-Myb to the KIX Domain of CREB

The KIX domain of CBP is a transcriptional coactivator. Concomitant binding to the activation domain of proto-oncogene protein c-Myb and the transactivation domain of the trithorax group protein mixed lineage leukemia (MLL) transcription factor lead to the biologically active ternary MLL∶KIX∶c-Myb complex which plays a role in Pol II-mediated transcription. The binding of the activation domain of MLL to KIX enhances c-Myb binding. Here we carried out molecular dynamics (MD) simulations for the MLL∶KIX∶c-Myb ternary complex, its binary components and KIX with the goal of providing a mechanistic explanation for the experimental observations. The dynamic behavior revealed that the MLL binding site is allosterically coupled to the c-Myb binding site. MLL binding redistributes the conformational ensemble of KIX, leading to higher populations of states which favor c-Myb binding. The key element in the allosteric communication pathways is the KIX loop, which acts as a control mechanism to enhance subsequent binding events. We tested this conclusion by in silico mutations of loop residues in the KIX∶MLL complex and by comparing wild type and mutant dynamics through MD simulations. The loop assumed MLL binding conformation similar to that observed in the KIX∶c-Myb state which disfavors the allosteric network. The coupling with c-Myb binding site faded, abolishing the positive cooperativity observed in the presence of MLL. Our major conclusion is that by eliciting a loop-mediated allosteric switch between the different states following the binding events, transcriptional activation can be regulated. The KIX system presents an example how nature makes use of conformational control in higher level regulation of transcriptional activity and thus cellular events

Detailed Regulatory Mechanism of the Interaction between ZO-1 PDZ2 and Connexin43 Revealed by MD Simulations

The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation